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Abstract 

We study the entanglement dynamics of strongly driven atoms off-resonantly coupled with cavity fields. 
We consider conditions characterized not only by the atom-field coupling but also by the atom-field de- 
tuning. By studying two different models within the framework of cavity QED, we show that the so-called 
atomic entanglement sudden death (ESD) always occurs if the atom-field coupling larger than the atom-field 
detuning, and is independent of the type of initial atomic state. 
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1. Introduction 

Entanglement plays a crucial role in quantum information processing [1]. Quantum algorithms 
(particularly, in Shor's algorithm, to find the prime factors of an n-bit integer) exploit entangle- 
ment to speed up computation. Entanglement dynamics have been a difficult subject and have 
attracted extensive interest recently ranging from two-qubit systems [2 — 5], to continuous vari- 
ables [6,7], spin systems [8,9,10], and multi-partite systems [11 — 14]. Entanglement of open 
quantum systems also has attracted considerable attention due to its significance for both funda- 
mentals and applications of quantum information processing [15, 16]. Such as, the paper [15] gives 
us the non-Markovian entanglement evolution of two uncoupled qubits which are coupled to their 
own independent environments. Origination and survival of entanglement for two independent 
qudits dephasingly coupled to a common zero-temperature, super-Ohmic, bosonic environment is 
studied in [16]. Moreover, in the close systems, proposals [17, 18] have been made for the direct 
measurement of finite-time disentanglement in cavity QED, and real-time detection of entangle- 
ment sudden death (ESD) has been reported very recently[19, 20, 21]. 

The interaction of a two-level atom with a quantized single mode of a harmonic oscillator, was 
called the Jaynes-Cummings (JC) model [22] . The JC model has found its natural playground in 
the field of cavity quantum electrodynamics (CQED), and extensions of the JC model to more 
atoms and more modes, externally driven or not, have been developed. Presently, we enjoy a vast 
number of theoretical and experimental developments. In [19], the authors dealt with a double 
JC model in which two initially entangled two-level atoms A and B are independently coupled 
with separate cavity fields a and b respectively, but there are no interactions at all between the 
subsystems Aa and Bb. Focusing on the atomic subsystem they found that depending on the type 
of initial state of atoms AB, their entanglement may or may not exhibit ESD. From [19], we can 
acquire that the authors have taken Bell-like states as the initial state of atoms and assumed the 
initial cavity field to be in vacuum. As a consequence, ESD was found to be sensitive to the initial 
atomic state. That is to say, ESD may occur for a certain type of initial atomic state but does not 
appear for another type. 

In this paper, we consider a system consisting of a two-level atom trapped inside a single 
mode cavity. The atom is additionally driven by a strongly classical field. The experimen- 
tal implementation seems to be feasible due to the recent advances in deterministic trapping of 
atoms in the optical cavities [23,24]. We show that under-taking atomic pure Bell-like states 
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|$ab) = cos6\e A ,e B ) + sin Oe^lgA, 9b) or \V A b) = cos6\e A ,g B ) + sin0e^|#A, e B ) (with 
< 9 < 2n and < < n) as the initial state of atoms and assuming the initial cavity field 
to be in vacuum, ESD can occur for two types of initial atomic states. In order to study entan- 
glement we will use the negativity (N), which can be defined for two qubits as two times the 
modulus of the negative eigenvalue of the partial transposition of the state p, o Ta [25, 26], if it 
exists. For short 

N{p) = 2 max{0, -A min }, (1) 

where X min is the lowest eigenvalue of o Ta . Our choice is motivated by the facts that the neg- 
ativity is easy to calculate and provides full entanglement information for a two-qubit system. 
We also need to quantify three-party entanglement in the tree-qubit state \ipABa), the three-tangle 
t^abo) has been introduced in Ref.[27]. It can be expressed by using the wave function coeffi- 
cients {^ooo, V'ooi, • • • , } as 

r 3 = A\d 1 -2d 2 + 4:d 3 \, 

d i = V'oooV'm + V'ooiV'iio + ^010^101 + ^ioo^on> 

d 2 = V'ooo^iii^oii^ioo + ^ooo^mV'ioi^oio + ^ooo^m^no^ooi 

+ V'Oll^lOO^lOl^OlO + ^011^100^110^001 + ^101^010^110^001, 

4 = V'oooV'iioV'ioiV'oii + V'ni^ooiV'oio^ioo- (2) 

To figure out general conditions for the possible negativity of atomic ESD, we shall study two 
different atom-cavity models, which we refer to as Model 1 and Model 2. Model 1 is considered 
in section 2, in this model one of two atoms is trapped in a single cavity, off -resonantly coupled to 
this cavity, and driven by a classic strong coherent field, while the other remains outside the cavity 
and has no environment. In section 3 we consider Model 2 which deals with a double driven JCM: 
each of the two strongly driven atoms interacts with its own cavity in the absence of any coupling 
between the atom-field subsystems. A generic result we have found out is that the atomic ESD 
always occurs in the conditions of the atom-field coupling larger than the atom-field detuning. 
Finally, we conclude in section 4. 

2. Model 1 

In this section we show a model in which a strongly driven atom A is off -resonantly coupled 
to a single-mode cavity field a, while the other atom B is isolated from all environment. The 
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Hamiltonian of the system can be described by 

H = ^y^A + ^y°b + ^/« f a + nO(e-^V+ + e iuJDt a A ) + Kg(o\a + a A a)), (3) 

where Vt is the Rabi frequency associated with the coherent driving field amplitude, g is the atom- 
cavity mode coupling constant, a{a)) the field annihilation (creation) operator, a A = \qa){^a\ 
{g\ = \e A )(gA\) the atomic lowering (raising) operator, and a A = \e A )(e A \ — \9a){9a\ the inver- 
sion operator. Considering the strong-driving regime for the interaction between the atom and the 
external coherent field Q^>{g, 5} and making uj a = uj d , we can use the rotating-wave approxima- 
tion (RWA) obtaining the effective Hamiltonian [28] 

Heff = y K + *A)(ae m + Je- m ). (4) 

where 8 = Uf — u a is the atom-cavity detuning. 

We first choose the state \^ab) as an initial state of atoms A and B, and assume the initial 
cavity field to be in vacuum (\Vacuum) means no photon in the single-mode cavity field). Thus 
the total system state at t = is 

|*(0)}ab o = (cos 9\e A ,g B ) + sin9e t(t, \g A ,e B )) <g> \Vacuum) a . (5) 

For simplicity, we make = 0. Then the evolved state in time t will be \^?(t)) = U (t)\^f(0)), and 
U(t) = exp(—iH e fft), so 

. T . . . cos 9 + sin 9 . . . . cos 9 — sin 9 . . . . 

\*(t))ABa = 2 \<*)\+A,+b) + 2 l«)l+A,-fl) 

cos 9 — sin 9 cos 9 + sin 9 
= \-a)\-A,+ B ) ^ \~ ot)\-A,-B), (6) 

with a = §(1 - e iSt ) and \± x ) = (\g x ) ± \e x ))/V2, where X = A, B and {\i)}j =1 = {\ + A 
+b), I +a ~b), I — a +b), I — a — b)} is the rotated basis of the atomic Hilbert space. We now 
try to estimate the atomic entanglement of the state Equation©, and we notice that the cavity 
field in the general nonorthogonal coherent state \a) and | — a). So we can define |0) = \a), 

2 

|1) = (| — a) — P\a))/yl — P 2 with P = exp(^(cos5t — 1)) for the cavity field subsystem. 

Then Equation© can be written as 

■ * ii . cos 9 + sin 9 . . cos 9 — sin 9 . 
W)> = ( 2 \+^+b) + 2 

cos 9 — sin 9 cos 9 + sin 9 
- 2 I ~ A +B ' 2 ~ A _B 

. cos 9 — sin r- — . . cos 9 + sin 9 



vt^i a + B ) + I v r^i -A -s»|l> (V) 







FIG. 1: N(pJ LB (t)) as a function of 6 and the dimensionless time gt for g/5 = 2, with 5 = 0.5g in Model 
1. 

In order to calculate the negativity of atoms, we trace out the mode a, the expression for p\ B (t) in 
the rotated basis of the atomic Hilbert space 

,cos0 + sin0, 2 ,cos9 — sin^, 2 

Pn = p 4 4 = ( ^ ) > p22 = P33 = ( 2 ' ' 

cos 9 + sin 9 2 cos^ — sin^ 2 
Pu = P4i = -( g ) P ' P23 = P32 = ~( 2 ^ 

cos 2 - sin 2 6> 

Pl2 — P21 — P34 = P43 — ^ j 

cos 2 9 — sin 2 9 

Pn = P31 = pA2 = P24 = -. P- (8) 



The negativity of p\ B {t) is 



N(p\ B {t)) = 2max{0, - P - y/l + P 2 - 2Pcos4^)}. 



(9) 



The negativity N(p\ B (t)) has been calculated, as shown in Fig.l and Fig. 2, for the different 
g/5. Obviously the negativity is fluctuating with dimensionless time gt and the mixing angle 9, and 
it can be zero in a finite amount of time (the so-called entanglement sudden death). An important 
point is that ESD is sensitive to the value of g/5 as displayed in Fig.l and Fig. 2. For g/5 > 1, 




FIG. 2: N(p^ B (t)) as a function of 6 and the dimensionless time gt for g/5 = 0.5, with 5 = 2g in Model 
1. 

ESD happens readily (Fig.l), and the larger the value of g/5 is, the more easily the ESD appears. 
When g/5 < 1, ESD will not appear at any time (Fig. 2). This phenomenon shows that ESD is 
related to g/5 completely in this model. Moreover, As we study the entanglement dynamics under 
a close system, Fig.l also shows that sudden death and resurrection of atomic entanglement can 
occur periodically as well as alternately. After entanglement of two atoms occurs sudden death for 
a finite time, atomic entanglement can revive again gradually. We have verified that when the ratio 
g/5 is larger than the one used in Fig.l, the region of ESD and entanglement revivals increases, 
and they last for a longer time. On the other hand, when we decrease g/5, the ESD regions shrink 
and the period of the entanglement oscillations becomes smaller. 

The character revealed by Fig. 3 is that 5 influences the period of the ESD and g is related to 
the velocity of the AB subsystem's disentanglement. If we want to acquire a long-time ESD, then 
we can choose a smaller atom-field detuning <5(Fig.3(a)). While 5 is a constant, the AB subsystem 
can rapidly disentangle with the larger g (Fig. 3(b)). A physical interpretation of the result is that 
the more the atom to the cavity mode couples, the easier the initial entanglement decays. Then we 
consider the atom and the cavity under resonance condition (5 = 0) or the large detuning regime 
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FIG. 3: (a) N(p%(t)) as a function of time t for g = 4, 9 = vr/4 in Model 1, (b) N(p\ B (t)) as a function 
of time t for <5 = 1, = 7r/4 in Model 1. 

(5^>g). Firstly, in the resonant regime the negativity of p\ B {t) decreases at the beginning, and 
then vanishes at any time and resurrection of atomic entanglement can not happen forever under 
this condition(Fig.4), so the strongly driven model under the resonant regime can be considered 
as an information eraser. Secondly, entangled states are generally very fragile against interaction 
with environments, but in our present scheme, the energy exchange between atom and cavity do 
not exist under the atom-cavity large detuning regime, so the initial entanglement of atoms can be 
preserved during system evolution, which is shown as Fig. 5. That is to say, the strongly driven 
and large detuning model can be used to preserve the information (entanglement) in quantum 
processing and quantum computing. 
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FIG. 4: N(p^ B (t)) as a function of 6 and the dimensionless time gt for 5 = in Model 1 

Now we move to the negativity iV(p* a (t)) and N(p Ba (t)) and the three-tangle T~{p\ Ba {t)), and 
we acquire the following 

iV(pI(t)) = 2max{0,Q^}, (10) 

iV(pf a (t)) = 2max{0,QBa}, (11) 

r(pi Ba (t)) = ±(l-P 2 )(l-cos46), (12) 

if < cos2# < 1, then Q Aa = -|(2 - 2cos2# - v 7 ^ 3 + (4 - 8P 2 )cos2# + cos4#) 
and Q Ba = -|(1 - Pcos2fl), and if -1 < cos2fl < 0, then Q Aa = -|(2 + 2cos2^ - 
V2y/3 + (8P 2 - 4)cos20 + cos4#) and Q Ba = -|(1 + Pcos29). Due to < P < 1, then 
QBa < 0, the negativity N(p Ba (t)) = 0. In the strongly driven regime, B is isolated from all 
environment, and the strong classical field that drives the atom A can be used to enhance the 
atom( A) -cavity field(a) interaction. The entanglement between A and cavity a can appear during 
the system evolution, but B and cavity a can not entangle with each other in the strongly driven 
regime(that is different from the standard J-C Model [21]). It is interesting to find that while the 
AB subsystem appears the long-time ESD, the entanglement of the Aa subsystem and the three- 
tangle of the whole system are all on the long-time invariable entanglement(as shown in Fig. 6(a)). 
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10° 



FIG. 5: N(p AB (t)) as a function of 9 and the dimensionless time gt for 5 = lOg, in the large detuning 
regime in Model 1. 

It is clearly shown that the three-tangle r ABa can be influenced by the initial state |^ab(0)) in 
Fig. 6(6). On the condition 9 = n /A (the initial state is on the maximum entanglement), when 
AB subsystem's entanglement vanishes for a long time, the three-tangle of the whole system can 
achieve the maximum value 1 at the same time and the entanglement between A and a can not 
appear for any time. That is to say, the AB subsystem's entanglement is transferred to the whole 
system's entanglement thoroughly (Fig. 6(c)). So, we can acquire the three-partite system's long- 
time maximum entanglement in the strongly driven regime (In Model 1). 

Next, we choose \®ab) to be the initial atomic state. Then the initial total system state is 

$ABa(0) = {cos6\e A , e B ) +sin0e'*|&4,0B)) ® \ Vacuum) a , (13) 
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with = 0, and for p\ Ba (t) we have the expressions as follows 

„ cos 9 + sin 9 2 cos $ — s i n & i 

),P22 = P33 = ( ^ 

9 ^ / cos — sin 9 , g „ 

) 2 P,P23 = P32 = ( ^ ) 2 P, 



Pll 


= P44 


= l 2 






cos 9 + sin 9 


Pl4 


= P41 


" ^ 2 


Pl2 


= P21 


— P34 — P43 — 


Pl3 


= P31 


— P42 — P24 — 



cos 2 9 - sin 2 9 
4 ' 
cos 2 9 — sin 2 9 

: P- (14) 
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The negativity of p AB (t) is satisfied with 

N(PA B (t)) = N(p* AB (t)). (15) 

Equation(fT5l) reveals that the behavior of N(p AB (t)) dependent on g/5 is equal to that of 
N(p AB (t)). We also calculate the negativity N(p Aa (t)) and N(p Ba (t)) and the three-tangle 

r (P*Ba W)> and find that the y a« a11 e q ual t0 ^(p^t)), N(p% a (t)) and r(p% a (t)). 

Consequently, for Model 1, where only one driven atom interacts with its environment, if the 
atom is not driven and the cavity field is in the vacuum initially, ESD does not occur either for 
Pab(O) or Pab(O)- m contrast, however, if the atom is driven by a strong classical field and both 
the atom-cavity coupling and the atom-cavity detuning satisfy g/5 > 1, the atomic subsystem AB 
always evolves via ESD, independent of the type of initial atomic state(i.e. it never matters if it is 
PSfl(O) orp^(O)). 

3. Model 2 

The double JCM, which has been intensively investigated recently [19], has shown that if atoms 
A and B are prepared in the |<&AB)-type Bell-like pure state initially, then the atomic entanglement 
dies in a finite time and remains dead for some time before reviving itself again, i.e., ESD occurs, 
whereas the \^ ab)-^V q Bell-like pure state dose not exhibit ESD at all. However, these investi- 
gations are confined to the atom-field under full resonance conditions. Our purpose here is focus 
on the differences that may appear on the conditions the atoms A and B independently driven by 
a strongly external classical field and the atom-cavity off resonance. 

In Model 2, we consider two remote two-level atoms A and B which are first prepared to be in 
an entangled state, and then let each atom couple with a single-mode cavity respectively. During 
the interaction with the single-mode cavity, the two atoms are independently driven by a strongly 
classical field. In the strong-driving regime, The effective Hamiltonian governing Model 2 is of 
the form 

R eff = E H «f> < 16 > 
k=A,B 

where 

H eff = ^rK + + *k){a k e«* + a\e-^). (17) 
The total system state at t = of the form 

\^(0)) ABab = (cos9\e A ,g B ) + sin9e l4> \g A ,e B )) <g> \ Vacuum) a \Vacuum) b . (18) 
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also with = 0, then the evolved state in time t will be 



\V{t))ABab = exp{-iH eff t)\V{0)) A Bab 

cos 9 + sin 9 cos 9 — sin 9 . 

= g \ a a,Pb)\+A, +B) H g \<Xa,—Pb)\+A,-B) 

cos6>-sin#. . . cos9 + sin9. . 

= 1 - +b) = 1 - "a, -Pfc)|-A, -B), (19) 

where {|i)}| =1 = {| +a+b},\ +a —b), I -a +b), I -a — b)} is the rotated basis of the atomic 
Hilbert space. Similar to Equation©, we also define 



|0 a > = \aa), |l a ) = (| - a a ) - P A \a a ))/Jl - Pi, (20) 



|0 6 ) = | - A), |1 6 > = -P B \- Pb))/y/l-Ph (21) 
where Pa = exp (§# (cos t^t— 1)),Pb = exp(f#Ycos <5et — 1)). The reduced density matrix p? B (£) 

A °B 

is 

. cos 9 + sin 2 cos 9 — sin 6* 2 

pll = p 4 4 = ( ~ ) ,P22 = P33 = ( ) , 

cos6» + sin6) 2 cos 9 - sin 9 2 

Pl4 = P41 = -( - ) PaPb,P23 = P32 = ~( ~ ) PaPb, 

cos 2 9 - sin 2 9 
P12 = P21 = P34 = P43 = ^ Pb, 

cos 2 g - sin 2 9 

P13 = P31 = P42 = P24 = ^ Pa- (22) 

The negativity of p\ B (t) is 

N(As(t)) = 2max{0, -\{2-2P A P B -y/2J '(1 + P 2 )(l + P 2 ) + (P 2 + P 2 - 1 - AP A P B - P 2 P 2 ) cos40)}. 



(23) 

The negativity of p* B (i) is given by Equation(l23l) in full similarity with Model 1, ESD always 
occurs for g A /5 A > 1 and gs/^B > 1- The physical interpretation of the result is that g A /5 A > 1 
(or gs/^B > 1) means the strong coupling regime. Because the atom A (or B) and cavity a 
(or b) couple with each other strongly, the entanglement of the AB subsystem can occur ESD 
phenomenon. g A and g B are related to the velocity of the initial state's disentanglement, and 5 A , 
5 B influence the ESD's period. In addition, the ESD also appears on the condition g A /5 A > 1, 
5b>9b or g B /8 B > 1, S A >g A , even when 5 A > g A , 5 B > g B (but when 5 A > g A or 5 B > g B , 5 A 
and 5b do not correspond to the large detuning regime), and is shown as Fig. 7 and Fig. 8, in which 
the ESD phenomenon can also occur in the case when 5 A and 5 B are larger than g A and g B (but not 
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FIG. 7: N(p^ B (t)) as a function of 9 and the dimensionless time gt for 5a = O.bg, 5b = 2g in Model 2. 




FIG. 8: N^p^gft)) as a function of 6 and the dimensionless time gt for ^ = <5b = !•% in Model 2. 
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by too much). As we all know, the larger atom-cavity detuning means that the interaction between 
the atom and cavity becomes small, for example, in the large detuning regime, the cavity and the 
atom can not exist energy exchange. So when 8a and 5 B are much larger than g A and g B , the 
ESD will disappear. If we consider the driven atoms and the cavities under resonance condition 
(8 A — 8 B — 0), the negativity of p AB (t) decreases at the beginning, and then vanishes for all time 
(the same as Fig. 4). So Model 2 under the resonance regime can also be used as an information 
eraser. 

For another type of atomic initial state \<& AB ), it can be verified that the negativity N(p\ B (t)) 
is also equal to N(p\ B (t)). Through calculating the negativity N Aa , N Bb , N Ab and N Ba in Model 
2, we find that the entanglement of Aa and Bb can occur during the system evolution and can be 
on a long-time invariable entanglement when the AB subsystem is occurring the ESD. Because of 
no existing interaction in the Ba subsystem (or Ab subsystem), it is impossible to entangle Ab (or 
Ba) in the strongly driven regime (these results are not similar to the standard double J-C Model 
[19]). Hence, our conclusion regarding Model 2 is that the driven atomic subsystem AB always 
suffers ESD if the atom-cavity coupling and the atom-detuning are satisfied with the condition 
g A /5 A > 1, g B /5 B > 1 or g A /5 A >l, 5 B >g B or g B /5 B >l, 5 A >g A , even when 8 A > gA, 8 B > g B 
(but 5 A and 5 B can not correspond to the large detuning regime), independent of the type of the 
atomic initial state which may be either p AB (0) or p* B (0). 

4. Conclusion 

We have described the entanglement evolution of two two-level atoms off-resonantly coupled 
to cavity fields. In Model 1, one of two atoms is trapped in a single cavity, off-resonantly cou- 
pled to this cavity, and driven by a classic strong coherent field, while the other remains outside 
the cavity and has no environment. However, in Model 2, each of the two strongly driven atoms 
interacts with its own cavity in the absence of any coupling between the atom-field subsystems. 
There are different available or forthcoming routes to the implementation of our model. In the 
microwave regime of cavity QED, pairs of atoms excited to Rydberg levels cross a high-Q super- 
conductive cavity with negligible spontaneous emission during the interaction [29]. In the optical 
regime the application of cooling and trapping techniques in cavity QED [23] allows the determin- 
istic loading of single atoms in a high-finesse cavity, with accurate position control and trapping 
times of many seconds [24] . In this regime laser-assisted three-level atoms can behave as effective 
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two-level atoms [30] . On the other hand, trapped atomic ions can remain in an optical cavity for an 
indefinite time in a fixed position, where they can couple to a single mode without coupling rate 
fluctuations [31]. These systems are quite promising for our purposes and could become almost 
ideal in the case of achievement of the strong coupling regime. 

Under off -resonance conditions and starting from the vacuum state of the cavity fields, for neg- 
ligible atomic decays and cavity leakage, we solved exactly the system dynamics for two types of 
initial preparation of the atom pairs (|<E>ab) and \ ^ab))- Thus we found conditions for the negativ- 
ity of the so-called atomic ESD. Namely, the initial entanglement of atoms, if any, will eventually 
suffer a sudden death, if the atom-cavity coupling and the atom-cavity detuning are satisfied with 
g/5 > 1 (in Model 1), while in Model 2 the system is accordance to Qa/^a > h 9bI$b > 1 
or 9a/$a > I, S B > g B and/or Qb/^b > 1, 5a > Qa (but 5 A and 5 B can not correspond to the 
large detuning regime). Furthermore, it is interesting to note that such conditions for ESD do not 
depend on the concrete type of the initial state of atoms. The atom-cavity detuning 5 influences 
the period of ESD and the atom-cavity coupling g is related to the velocity of the AB subsystem's 
disentanglement. 

In the atom-cavity resonance regime, the negativity of the atomic subsystem decreases at the 
beginning, and then vanishes for any time, so the strongly driven model (such as Model 1 and 
Model 2) under the resonance regime can be used as an information eraser. Entangled states are 
generally very fragile against interaction with environments, but in our present scheme, the energy 
exchange between atoms and cavities do not exist under the atom-cavity large detuning regime, so 
the initial entanglement of atoms is preserved during system evolution. That is to say, the strongly 
driven and large detuning model can be used to preserve the information (entanglement) in quan- 
tum information processing and quantum computing. 

In the strongly driven regime, we note the negativity of two atoms is a periodic function of 
time in the close system. It means that atoms evolve quasiperiodically between the entangled and 
disentangled states. This result is different from refs.[15, 16], ref.[15] shows that the negativity of 
the bipartite system is a mononically decreasing function of time when the first qubit is trapped 
in a infinite heat bath and the second qubit has no environment. But authors in ref.[15] also have 
proposed that in the case when the second qubit is set in a finite and controlling quantum environ- 
ment, the dynamics of entanglement can oscillate, this nonmontonic behavior appears to be linked 
to the non-markovian character of the dynamics. That is to say, non-markovian entanglement can 
evolve quasiperiodically between the entangled and disentangled states in the open systems. From 
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the paper [15] we acquire the degree of entanglement at certain time intervals can be manipulated 
by means of an appropriate choice of the initial state of the controlling quantum environment. Our 
results in the strongly driven regime manifest the entanglement of two atoms can be manipulated 
through controlling the atom-cavity detuning and atom-cavity coupling. 
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